Numerical Stability Analysis of an Acceleration Scheme for Step Size Constrained Time Integrators
نویسندگان
چکیده
Time integration schemes with a fixed time step, much smaller than the dominant slow time scales of the dynamics of the system, arise in the context of stiff ordinary differential equations or in multiscale computations, where a microscopic time stepper is used to compute macroscopic behaviour. We discuss a method to accelerate such a time integrator by using extrapolation. This method extends the scheme developped by Sommeijer (Comput. Math. Appl. 19 (6):37–49 1990), and uses similar ideas as the projective integration method. We analyse the stability properties of the method, and we illustrate its performance for a convection-diffusion problem.
منابع مشابه
On the Construction of Variational Integrators for Optimal Control of Nonholonomic Mechanical Systems
In this paper we derive variational integrators for optimal control problems of nonholonomic mechanical systems. We rewrite the system as a constrained second-order variational problem, that is, as a problem where the Lagrangian and constraints are defined in terms of the position, velocity and the acceleration of the system. Instead of discretizing directly the equations of motion, we discreti...
متن کاملNonlinear Stability Analysis of Area-Preserving Integrators
Linear stability analysis is inadequate for integrators designed for nondissipative systems such as Hamiltonian systems in which nonlinear effects are often decisive. Mathematical theory exists (KAM theory) for rigorous analysis of small perturbations from equilibria, but it needs to be expressed in a form that is more easily applicable to the study of area-preserving maps. We have pursued this...
متن کاملStability of Asynchronous Variational Integrators
The formulation of multiple–time–step integrators can provide substantial computational savings for mechanical systems with multiple time scales. However, the scope of these savings may be severely limited by the range of allowable time step choices. In this paper we have performed an exhaustive study of the linear stability of the fully asynchronous methods called AVI (asynchronous variational...
متن کاملHigh-Order Symplectic FDTD Scheme for Solving Time-Dependent Schrödinger Equation
Using the three-order symplectic integrators and fourth-order collocated spatial differences, a high-order symplectic finite-difference time-domain (SFDTD) scheme is proposed to solve the time-dependent Schrödinger equation. First, the high-order symplectic framework for discretizing Schrödinger equation is described. Then the numerical stability and dispersion analyses are provided for the FDT...
متن کاملChallenges in Geometric Numerical Integration
Geometric Numerical Integration is a subfield of the numerical treatment of differential equations. It deals with the design and analysis of algorithms that preserve the structure of the analytic flow. The present review discusses numerical integrators, which nearly preserve the energy of Hamiltonian systems over long times. Backward error analysis gives important insight in the situation, wher...
متن کامل